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We present a microscopic theory for the d.c. Josephson ef- 
fect in model mesoscopic constrictions. Our method is based 
on a non-equilibrium Green function formalism which allows 
for a self-consistent determination of the order parameter pro- 
file along the constriction. The various regimes defined by the 
different length scales (Fermi wavelength Af, coherence length 
£o and constriction length Lc) can be analyzed, including the 
case where all these lengths are comparable. For the case 
\f<(Lc,£o) phase oscillations with spatial period Af/2 can 
be observed. In the case of Lc > £o solutions with a phase- 
slip center inside the constriction can be found, in agreement 
with previous phenomenological theories. 



PACS numbers: 74.50.+r, 85.25.Cp, 73.20.Dx 

The problem of d.c. Josephson transport through a su- 
perconducting constriction received a considerable atten- 
tion several years ago in the context of what was known 
as the theory of superconducting weak- links [j]. In the 
last few years, there has been a renewed interest on the 
subject fostered by the advances in nano-scale technolo- 
gies which would allow the fabrication of superconduct- 
ing mesoscopic devices. Actually, some steps in this di- 
rection have already been taken with the development of 
a Josephson Field Effect Transistor (JOFET) §. An- 
other experimental situtation where this problem can be 
studied is that of a superconducting Scanning Tunneling 
Microscope ||. 

The relevant feature of such a mesoscopic device would 
be the phase coherence of both single-electrons and 
Cooper pairs over a length comparable to the system 
size. This opens the possibility of observing novel in- 
terference phenomena when < (£o> £jv> Lc) (£jv de- 
notes the normal electrons coherence length) . The case in 
which Lc <C £o has been analyzed in a recent publication 
by Beenakker and van Houten [Q . In this regime the de- 
tailed form of the order parameter profile inside the con- 
striction is irrelevant for the evaluation of the Josephson 
current. In the opposite limit (Lc ^ £o) the supercon- 
ducting phase would drop linearly along the constriction. 
A model for a superconducting point contact within this 
limit has been proposed by Furusaki et al ||. 

In the intermediate regime, in which all the rele- 
vant lengths can be comparable, one would need a self- 
consistent determination of the complex order parameter 
along the constriction. This poses a difficult task that has 
only been addressed within phenomenological Ginzburg- 
Landau theory M or with semi-classical treatments based 



on Boltzman-type transport equations j7]]. The aim of 
this work is to present a microscopic model for the self- 
consistent description of a Josephson constriction in all 
the relevant regimes. 

We consider a model constriction like the one depicted 
in Figure 1, which consists on a quasi one dimensional 
region of length Lc coupled to wider regions (L and 
R) which act as electron reservoirs. For simplicity, we 
assume a single quantum channel open for transport 
through the constriction (a generalization for the multi- 
channel case is possible within the formalism outlined 
below) . The reservoirs are homogeneous superconductors 
with a constant complex order parameter except near the 
interface region. Within our model, the constriction may 
be either superconducting or normal, although we will 
focus here on the superconducting case. 

We use a local orbital representation for the description 
of the electronic states of the system. This representa- 
tion is very well suited for the self-consistent determina- 
tion of the order parameter profiles, provided that the 
pairing potential is assumed to be diagonal in this ba- 
sis. This simplification does not affect in a relevant way 
the description of the superconducting state, which is 
controlled basically by a single parameter (the coherence 
length Co ~ hvp/irA). In the lower part of Fig. 1 we rep- 
resent, schematically, our discretized model constriction 
where the quasi one dimensional region is represented by 
a linear chain with Nc sites. For convenience, in the 
numerical self-consistent calculations, the left and right 
reservoirs are represented by two Bethe lattices with co- 
ordination number z + 1. The mean- field hamiltonian for 
this model has the form 



(i) 



where the chemical potential fi is a constant throughout 
the whole system due to the absence of an applied voltage 
and will be taken as zero. The hopping parameters Uj 
are only different from zero for nearest neighbours and 
the complex order parameters are given by 



A, 



-Ui < 4^ >, 



(2) 



where Ui is the attractive e-e interaction at site i. The 
system defined by Eq. (2) provides the set of self- 
consistent conditions for obtaining the order parameter 



profile. We choose Ui and iy to be constant (Ul, Uc, Ur 
and tL, tc, £r) inside each different region (left reservoir, 
constriction, right reservoir) in such a way as to fix the 
desired bulk values of A (A L , A^, A° R ) in the three sepa- 
rate regions. The hopping parameters, Uj, can be chosen 
to be real and we take — t R — t as the unit of energy. 
We denote by t^c and tcR the parameters coupling the 
left and right electrodes to the constriction. 

In the absence of an applied bias voltage, the Joseph- 
son current through the constriction between sites i and 
i + 1 is given by 

/»(*) = J t C J2 (< 4r(*)Ci+l*(*) > - < 4fla(*)M*) > 

(3) 

One can show that is independent of the chosen site, 
i.e. the continuity equation is fulfilled, only when the 
solution for the different A^ is fully self-consistent |^|j9) . 

The averaged quantities appearing in Eqs. (2) and 
(3) are most conveniently expressed in terms of non- 
equilibrium Green functions G" 1 |jTo|| , which in a 
Nambu (2x2) representation arc defined by : 



< c 



Ut')cn(t) > < c n (t')c^{t) > 



< 



> 



(4) 



Then, the self-consistent Eqs. (2) and the stationary cur- 
rent, Eq.(3), are given, in terms of the Fourier-transform 



A, 



Ui_ 



/CO 
du[G±-(u)} 21 , 
-OO 



2e f°° 

h = T tc J du ([G£i>)]ii - [Gfr +1 (co)h 



(5) 



(0) 



All the different correlation functions appearing in 
Eqs. (5) and (6) can be obtained using conventional Green 
function techniques (jl2| . 

In the limit of Lc it is possible to obtain analyti- 
cal results for the Josephson current using this model. In 
this case, the self-consistent phase profile is well approxi- 
mated by a step function between the L and R electrodes. 
Defining <p — 4>l — 4>r as the total phase diference, Eq.(6) 
can be written as 



— t LR sin(0) 



dulm[ 



l/M- (?) 



where tLR is the coupling between the outermost sites 
on the left and right reservoirs, f(u>) is the Fermi dis- 
tribution function; g£ and g r R are the retarded Green 
functions for the uncoupled electrodes (the tilde indicates 



that the phase factor has been removed, i.e. g r L 21 = 
e^s£,2i and g R12 = e~^g RA2 ), and 

Z>» = det[I - tl R r 3S l(u;)T 3SR (u;)}, (8) 

T3 being the usual Pauli matrix. 

It is interesting to analyze the transition from the tun- 
nel to the contact regime as given by Eq. (7) (for sim- 
plicity we assume A L = A R — A, and A <C 1 )• In the 
tunnel regime, tLR <C 1, D r ~ 1, and one recovers the 
usual Josephson expression for a tunnel junction |13[ 



7rA(2"i) 

I(4>) = sin tanh 

2eitAr 



A(T) 



2k B T 



(9) 



where Rn is the normal resistance of the junction which 
is given by Rjj 1 = (2e 2 /h)a, a being the normal trans- 
mission through the constriction which in our model 
adopts the simple form a — {At 2 LR / zt 2 ) / (I + t LR /zt 2 ) 2 
pTj . On the other hand, when approaching the contact 
limit, a — ► 1, the main contribution to the current comes 
from states inside the superconducting gap, which are 
given by the zeros of Eq.(8). These states are originated 
by multiple reflection processes at the interface region 
and give the following simple expression for the current 



2eR 



N 



■sm(0)tanh[^-^] 



(10) 



e{4>) = ±|A(T)|Jl - asin 2 (</>/2) being the position of 
the states inside the gap . For the special case a = 1 Eq. 
(10) yields I(<j>) ~ sin(0/2), which coincides with the re- 
sult given in Refs. [[t|Jl5[|. It is remarkable that, when 
a — > 0, Eq.(10) tends exactly to Eq.(9), whose deduction 
involved no localized states. One can understand the 
equivalence between both ways of obtaining the tunnel 
limit by realizing that the localized states move towards 
the gap edges when a — > 0, gradually becoming the band 
edge singularities of the uncoupled system. These singu- 
larities gave the main contribution to the current when 
obtaining Eq.(9). 

Let us turn our attention to the effects of self- 
consistency as one moves from the above regime to the 
opposite case, Lc/£o ^ 1- We have performed calcula- 
tions for a wide range of constriction lengths and different 
values of the normal transmission parameter, a. The co- 
ordination number z on both electrodes has been taken 
equal to 3, which ensures a rapid convergence to the bulk 
values of the order parameter. The self-consistent order 
parameter profile is calculated in the following way: we 
first assume a given initial profile along the constriction. 
Then, we calculate the system Green functions for this 
profile which in turn yield the new order parameter pro- 
file using Eq. (5). This process is repeated until conver- 
gence is achieved. 

In order to clarify the Lc / £o dependence of our results 
we have fixed A° c = A° L = A° R = 0.05 and Xr = 4a, (a 



being the intersite distance), throughout the whole cal- 
culation. For the same reason we restrict ourselves to 
the zero temperature case in the present calculations. In 
Fig. 2 we show I versus <p curves chosen to represent typ- 
ical behaviours found in the different physical regimes. 
Thus, Fig. 2a corresponds to the case of maximum trans- 
mission (a = 1) , whereas Figs. 2b and 2c illustrate cases 
with decreasing a. In these figures we plot curves for 
different values of Lc/£o, ranging from Lc/£,o = to 
Lc/£,o 1- We only represent the part of the curves 
with I((f>) > 0. 

Several features are noticeable in these curves. First, 
when increasing Lc/^o, a critical value Lc/^o ~ 1 is 
reached above which the function I(<p) becomes multival- 
ued |jl|,|U9) . This situation corresponds to the appearance 
of a second kind, or "solitonic" , solution with a phase slip 
center inside the constriction, in addition to the normal 
one consisting basically on a linear phase profile between 
the reservoirs. In Figure 3 the two kind of self-consistent 
profiles for Lq > £o are displayed. They correspond to 
(j> ~ 7T where the phase drop inside the constriction is spe- 
cially sharp. It is worth noticing the appearance of oscil- 
lations of period A^/2 in both profiles (phase and modu- 
lus). This interference effect is a consequence of the phase 
coherence of the normal electrons along the constriction, 
analogous to the oscillations found in the electrochemi- 
cal potential in normal mesoscopic wires |l6|,|l7j ■ On the 
other hand, the other relevant length scale, £o> manifests 
itself in a clear way in the solitonic phase profile: one can 
verify that the width of the phase-slip center is almost 
equal to £o- The overall features of these solitonic solu- 
tions are in agreement with predictions made in previous 
phenomenological analysis JUJlfifl . 

It is interesting to observe the behaviour of the max- 
imum Josephson current, I max , as a function of Lq for 
the different situations represented in Fig. 2. In the 
case of maximum transmission (Fig. 2a) a slow decrease 
in Imax is observed for increasing Lc, with a limiting 
value corresponding to the critical current of the infinite 
one-dimensional chain. This monotonic decrease of I ma x 
is not found when the trasmission is lowered below a cer- 
tain value. For instance, for the case of Fig. 2b, I max 
slightly increases with Lc, tending to the same limiting 
value as in Fig. 2a. 

Fig. 2c corresponds to a somewhat different physical 
situation, in which a weaker coupling between the con- 
striction and the reservoirs is considered. For sufficiently 
large Lc, this situation may be regarded as a case not 
far from two identical tunnel junctions connected in se- 
ries, which would have a I(4>) ~ sin(0/2) characteristic, 
whereas for Lc < £o the system behaves like a single 
junction with the usual I{4>) ~ sin0 form. The transi- 
tion between these two regimes can be clearly observed 
in fig. 2c. Note that with increasing Lc, I{(j>) becomes 
eventually multivalued, as in the previous cases. How- 
ever, for large enough Lc (see case Nc — 80 in fig. 2c), 



the solitonic branch does not merge in a continuous way 
into the "normal" one. Instead, both branches extend 
up to = 27r keeping a different character. In the lower 
branch, the region where the modulus of the order param- 
eter nearly goes to zero tends to fill the whole constriction 
when 4> — > 2tt. 

In conclusion, we have presented a fully self-consistent 
description of the d.c. Josephson transport on a meso- 
scopic weak-link. As a first check, the method has been 
applied to the case of a short constriction, recovering pre- 
viously known results. On the other hand, we have ana- 
lyzed the case where the constriction length is compara- 
ble or larger than the superconducting coherence length. 
To our knowledge, this is the first calculation of the self- 
consistent order parameter profiles and the current-phase 
relationship based on a microscopic model. The local 
character of our procedure makes it specially powerful 
for treating systems with a more complicated geometry 
in situations where the self-consistent variation of the 
order parameter plays a crucial role. For instance, us- 
ing the present approach one would be able to address 
the problem of the proximity induced Josephson effect in 
a metal-superconductor nano-junction |l9| . Work along 
this line as well as on the inclusion of a finite voltage for 
the study of the a.c. Josephson effect is under progress. 
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FIG. 1. Schematic representation of the model constriction 
considered in this paper. 

FIG. 2. Josephson current versus total phase difference for 
three different values of the transmission parameter a: (a) 
a = 1, (b) a — 3/4, and (c) a ~ 1/12. In cases (a) and 
(b) the hopping parameter between the constriction and the 
reservoirs {the and ten) is equal to the intra-chain hopping 
tc {tc = \fz for case (a) and tc = 1 for case (b)), whereas 
in case (c) we take tc = 1 and the = tcR = 0.5. The num- 
bers above each curve indicate the number of sites within the 
constriction. 

FIG. 3. Self-consistent order parameter phase and modulus 
profiles for the case Nc = 64 and <j> = 3.4 of Fig. 2 (b). Solid 
and dotted lines correspond to the solution in the upper and 
lower branches of the I(<f>) curve respectively. 



